DISCRETE TOMOGRAPHY OF PLANAR MODEL SETS 
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Abstract. Discrete tomography is a well-established method to investigate finite point sets, 
in particular finite subsets of periodic systems. Here, we start to develop an efficient approach 
for the treatment of finite subsets of mathematical quasicrystals. To this end, the class of 
cyclotomic model sets is introduced, and the corresponding consistency, reconstruction and 
uniqueness problems of the discrete tomography of these sets are discussed. 



1. Introduction 

Discrete tomography is concerned with the inverse problem of retrieving information about 
some discrete object from (generally noisy) information about its incidences with certain query 
sets. A typical example is the reconstruction of a finite point set from its line sums in a small 
number m of directions. The term X-ray (or X-ray projection) is a generic name here which 
stands for a mechanism that produces weighted projection data. More precisely, a (discrete 
parallel) X-ray of a finite subset of Euclidean d-space IR"^ in direction u gives the number of 
points in the set on each line in M!^ parallel to u. (This concept should not be misunderstood 
in the sense of diffraction theory, where X-rays provide rather different information on the 
underlying structure that is based on statistical pair correlations; compare with Guinier (1994), 
Cowley (1995) and Fewster (2003).) 

Many papers focus on the discrete tomography of subsets of lattices since lattices are good 
models for crystalline structures. However, nature provides us also with structured non-lattice 
sets, the so-called quasicrystals. In the present paper, we shall investigate the discrete tomog- 
raphy of systems of aperiodic order, more precisely, of so-called model sets (or mathematical 
quasicrystals), which are commonly accepted as a mathematical model for perfect quasicrys- 
talline structures in nature (Steurer, 2004). As model sets possess a 'dimensional hierarchy', 
which means that any model set in d dimensions can be sliced into model sets of dimension 
d — 1, solving the reconstruction problem for two-dimensional systems with aperiodic order 
lies at the heart of solving the corresponding problem in three dimensions. 

The main motivation for our interest in the discrete tomography of model sets comes from 
the demand of materials science to reconstruct three-dimensional (quasi) crystals or planar lay- 
ers of them from their images obtained with quantitative high resolution transmission electron 
microscopy (HRTEM) in a small number of directions. 

In fact, in Schwander et al. (1993) and Kisielowski et al. (1995), the technique QUAN- 
TITEM (quantitative analysis of the information coming from transmission electron mi- 
croscopy) is described, which is based on HRTEM and can effectively measure the number of 
atoms lying on lines parallel to certain directions. At present, the measurement of the number 
of atoms lying on a line can only be achieved for some crystals; see Schwander et al. (1993) 

1 



2 



BAAKE, GRITZMANN, HUCK, LANGFELD, AND LORD 



and Kisielowski et al. (1995). However, it is reasonable to expect that future developments in 
technology will improve this situation. 

Roughly speaking, planar model sets are projections of certain subsets depending on some 
window of a higher dimensional lattice into the plane. In Section [31 we will define model 
sets in general, but we will mainly restrict ourselves to a well-known class of planar model 
sets, the cydotomic model sets. On the one hand, cyclotomic model sets exhibit a particularly 
nice and useful algebraic structure, while on the other hand real-world quasicrystals can be 
sliced into parallel planar layers that can be modeled by cyclotomic model sets (Pleasants, 
2000). Also, in a certain sense, cyclotomic model sets can be seen as a direct generalization 
of the square lattice , the classical planar setting of discrete tomography. 

Naturally, all classic issues of discrete tomography including uniqueness, reconstruction 
and stability (see e.g. the book by Herman & Kuba (1999) and, in particular, the papers by 
Gardner & Gritzmann (1997), (Gritzmann et al. (1998), Gardner et al. (1999), Gritzmann 
et al. (2000), Alpers et al. (2001) and Alpers & Gritzmann (2006)) can be studied for model 
sets as, in principle, they are just different ground sets, for the potential solutions. As it turns 
out, however, the more general setting does disclose some new aspects, and the present paper 
will stress these. In particular, it is a priori not even clear how to decide whether a translate 
of a given finite point set occurs within an aperiodic structure. 

As a matter of fact, previous studies have focussed on the 'anchored' case that the underlying 
ground set is located in a linear space, i.e., in a space with a specified location of the origin. 
The X-ray data is then taken with respect to this localization. This assumption is mainly 
justified by the fact that, as point sets, one has the equahty t + = foj. ^11 t £ Z^. 
Hence, in the lattice case, one can always assume that - if a solution exists - it is close to the 
origin. In the affine and aperiodic case of planar model sets it is a priori not clear how far out 
solutions may exist and how one can systematically search for them. 

The main result of this paper is, however, that for cyclotomic model sets (coming from 
polyhedral windows) all possible localizations can be determined efficiently. In fact, we shall 
solve a corresponding decomposition problem and a separation problem. This will allow us to 
reduce tomographic problems such as reconstruction and uniqueness for cyclotomic model sets 
to the corresponding classical problems with certain restrictions. One difference is manifest in 
the fact that potential solutions are subsets of a finite list of patches, whose number typically 
grows polynomially in the size. In fact, using the algebraic and the geometric structure of 
cyclotomic model sets we show that in a well-defined way the algorithmic methods that have 
been developed for the lattice case can be extended to the discrete tomography of cyclotomic 
model sets. (Note, however, as a warning that even in the (linear) lattice case these 
problems are NF-hard for three or more lattice directions; see Gritzmann et al. (1998) and 
Gardner et al. (1999).) 

Let us be more specific. By using the Minkowski representation of algebraic number fields, 
we introduce, for n ^ {1,2}, the corresponding class of cyclotomic model sets Tl C C = 
which live on C C, where Cn is a primitive nth root of unity in C, e.g., Cn = e^'^*/". 

(Here, and in the following a subset S of is said to live on a subgroup G of if its 
difference set S — S := {s — s' | s,s' G S} is a subset of G. Obviously, this is equivalent to 
the existence of a suitable t E such that S C t + G.) The Z-module Z[Cn] is the ring of 
integers in the nth cyclotomic field Q(Cn), and, for n ^ {1, 2, 3, 4, 6}, when viewed as a subset 
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of the plane, is dense; see Section [21 for details. In contrast, (cyclotomic) model sets A are 
Delone sets, i.e., they are uniformly discrete and relatively dense. In fact, model sets are even 
Meyer sets, meaning that also A — A is uniformly discrete; see (Moody, 2000). It turns out 
that, excepting the cyclotomic model sets living on Z[Cn] with n G {3,4,6} (these are exactly 
the translations of the square and the triangular lattice, respectively), cyclotomic model sets 
A are aperiodic, meaning that they have no translational symmetries. Well-known examples 
with A^-fold cyclic symmetry are the vertex sets of the square tiling {n = N = 4), the triangle 
tiling (2n = = 6), the Ammann-Beenker tiling {n = N = 8), the Tubingen triangle tiling 
(2n = N = 10) and the shield tiling {n = N = 12), respectively; see below for details. Observe 
that 5, 8, 10 and 12 are standard cyclic symmetries of genuine planar quasicrystals (Steurer, 



Whether or not one has future applications in materials science of quasicrystals in mind, the 
starting point will always be a specific structure model. This means that the specific type of the 
(quasi) crystal is known, and one is confronted with the X-ray data of an unknown finite subset 
of it. Let us point out that the rotational orientation of the probe in an electron microscope 
can rather easily be ascertained in the diffraction mode, prior to taking images in the high- 
resolution mode, though a natural choice of a translational origin is not possible. Hence a 
first task is to 'localize' a given probe within Z[(^„]. To be more specific, suppose X-rays of 
some planar (quasi) crystalline set F are taken in some directions oi,...,Om £ ^[Cn] \ {0}. 
Obviously, every point of F is 'registered' by every X-ray image, hence F is contained in the 
grid 



see Definition IHl Of course, in general G contains many more points than F, hence does not 
disclose F. On the other hand, only those subsets F' of G whose X-rays coincide with the 
given data are feasible solutions which lie in a translate of the underlying model set. Hence 
a first problem is to determine the decomposition of G into the subsets which are compatible 
with the underlying Z-module Z[Cn], i-e., which lie in a common translate of Z[Cn]; see Section 
for details. This problem has its origin in the practice of quantitative HRTEM since, in 
general, the X-ray information does not allow us to locate the underlying Z-module 
Using standard results of algebra, we will actually show much more, namely that the solution 
of this decomposition problem only depends on n and the given Z[C„]-directions but not on 
the specific X-ray data. Hence, conceptually, we can consider the different equivalence classes 
separately. 

Of course, even if a Z[(^„]-equivalence class of the grid G contains a set F' whose X-rays 
coincide with the given data, this set need not belong to the underlying cyclotomic model 
set. Hence it is clear that additional constraints that are induced by the construction rules of 
the underlying model set have to be satisfied to guarantee feasibility. One possible approach 
could be to first reconstruct a potential solution that is compatible with the given X-ray 
data and then check whether it actually belongs to the underlying model set. Unfortunately, 
this approach does not lead to an efficient algorithm (see Remark 1^ . Therefore, we use the 
specific structure of model sets (originating from some window through a projection process) 
and determine which subsets of G can possibly arise. In fact, all possible solutions (that might 
actually lie 'far out' in the defining model set) can be found and explored by translating the 
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given window; see Section 01 For many types of windows, this separation problem can be 
handled by geometric techniques based on the theory of arrangements; see Sectional 

The present paper is organized as follows. 

As a service to the reader, we begin with two preliminary sections that put together the 
notions required and recall several tools from algebra and the mathematical theory of qua- 
sicrystals. In fact, the algebra is needed not only to properly explain cyclotomic model sets 
but is crucial for devising algorithms for checking containment of points in this structure and 
also yields best-known bounds for the running time of our basic algorithms. (Of course, in 
view of the prominent role of group theory in crystallography and materials science, the rele- 
vance of algebraic methods for our cyclotomic structures does not really come as a surprise.) 
In Section [21 we explain the algebraic concepts in an elementary way while Section [31 gives a 
concise but sufficiently detailed account of model sets. In particular, we introduce the special 
class of cyclotomic model sets, which will be the central objects of the present paper. Some 
examples illustrate the structure and the beauty of cyclotomic model sets. 

The key problems and main results will be formulated in Section 0} their proofs will be 
given in Section 

2. Algebraic Background and Notation 

For all n G N, and Cn a fixed primitive nth root of unity in C {e.g., Cn = e^'^*/"), let Q(Cn) 
be the corresponding cyclotomic field, i.e., the smallest intermediate field of the field extension 
C/Q that contains Cn- Further, denoting by Cn the complex conjugate of Cn, it is well known 
that Q(Cn + Cn) (defined analogously) is the maximal real subfield of Q(Cn), i-e., 

Q(C„) n R = Q(Cn + Cn) ; 

see Washington (1997, p. 15). Throughout this text, we shall use the notation 

K„ = Q(Cn), k„ = Q(Cn + Cn), ©n = ^Cn], On = Z[C„ + (n] , 

where Z[Cn] (resp., Z[Cn + Cn]) is defined as the smallest subring of C that contains Z and Cn 
(resp., Z and Cn + Cn)- Further, cp will always denote Euler's phi-function (often also called 
Euler's totient function), i.e., 

(p{n) = card {{k € N | 1 < A; < n and gcd(/c, n) = 1}) . 

Occasionally, we shall identify C with R^. 

The set On hosts the corresponding cyclotomic model sets (c/. Section rT^ : K^i, kn, and 
On will be needed for the analysis of the algebraic structure of On that will allow the relevant 
algorithmic computations. The following lemma shows how On is related to o„. 

Lemma 1. For n> 3, one has: 

(a) On is an On-module of rank 2. More precisely, one has On = On + OnCn, o,nd {1, Cn} 
is an On-basis of On- 

(b) K„ is a hn-vector space of dimension 2. More precisely, one has K„ = k„ + kn Cn, 
and {l,Cn} is a hn-basis o/K„. 

Proof. First, we show (a). The linear independence of {1, Cn} over On is clear: by our assump- 
tion n > 3, {l,Cn} is even linearly independent over R. For the remainder of the assertion 
we prove that all non-negative integral powers Cn satisfy Cn ~ + /^Cn for suitable a, (3 &On- 
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Figure 1. Central patches of the square tihng (left) and triangular tiling (right). 



Using induction, it suffices to show Cn = + PCn for suitable a, /3 S On- To this end, note that 
Cn = Cn^ and observe that = -I + ((„ + C ^)Cn- 

Claim (b) follows similarly. □ 



Remark 1. Seen as a point set of R^, On has A^-fold cyclic symmetry, where 
(1) = iV(n) := lcm(n,2) = 



n, if n is even, 
2n, if n is odd. 

Except for the one-dimensional case n £ {1,2} (Oi = O2 = Z), the crystallographic cases 
n £ {3,6} (triangular lattice Os = Oq, see Figure 121 and n = 4 (square lattice ©4, see 
Figure 121), is dense in R^. For the latter, note that, by Lemma On is an o^-module of 
rank 2, whose R-span is all of R^. For n G N \ {1,2,3,4,6}, On is a Z-module of rank > 2 
(see Remark IHl below) embedded in R, hence a dense set in R. Consequently, On is then a 
dense set in R^. 

The following well-known result is needed later to actually compute the coordinates of 
On-points. As usual, denotes the group of units of a given ring R. 

Proposition 1 (GauiJ). One has [K„ : (Q] = (/)(n) and {1, Cn, Cn' • • • ' Cn*'"'' ^} is a Q,-basis 
o/K„. Moreover, the field extension K„/(Q is a Galois extension with Abelian Galois group 
G(K„/(Q) = {1i/n1i)^ , where a(modn) corresponds to the automorphism given by Cn ' — ^ Cn- 

Proof. See Theorem 2.5 of Washington (1997) and, for the statement about the Q-basis, the 
proof of Proposition 1.4 in Chapter V.l of Lang (1993). □ 

Remark 2. Note the identity 

{Z/nZy = {a(modn) | (a,n) = 1} 

and consult Table 3 of Baake & Grimm (2004) for examples of the explicit structure of 
G(K„/Q). 
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Corollary 1. If n > 3, one has [k„ : Q] = (j){n)/2. Moreover, a (Q,-basis o/k„ is given by the 

set {1, (Cn + Cn), (Cn + Cn? , • • • , (Cn + Cn)^^"^/'"' } • 

Proof. The statement about the degree [k„ : Q] is an immediate consequence of Lemma^b), 
Proposition n and the 'degree formula' for field extensions: If E/F/K is an extension of fields, 
one has [E : K] = [E : F][F : K] {cf Chapter V.l, Proposition 1.2 of Lang (1993)). The 
statement about the Q-basis again follows from the proof of Proposition 1.4 in Chapter V.l 
of Lang (1993). □ 

A full Z-module (i.e., a module of full rank) in an algebraic number field K which contains 
the number 1 and is a ring is called an order of K. It turns out that among the various 
orders of K there is one maximal order which contains all the other orders, namely the ring 
of integers in K; see Chapter 2, Section 2 of Borevich & Shafarevich (1966). For cyclotomic 
fields, one has the following well-known result. 

Proposition 2. For n e N, one has: 

(a) On is the ring of cyclotomic integers in and hence is its maximal order. 

(b) On is the ring of integers o/k„, and hence is its maximal order. 

Proof. See Theorem 2.6 and Proposition 2.16 of Washington (1997). □ 

Remark 3. It follows from Proposition ISl^a) and Proposition ^ that On is a Z-module of 

rank 0(n) with Z-basis {1, C„, C^, . . . , Cn^"^ ^}. Likewise, Proposition |2{b) and Corollary ^ 
imply that, for n > 3, is a Z-module of rank 4){n)/2 with Z-basis given by the set 

{1, (Cn + Cn), (Cn + Cn)^ • • • , (Cn + Cn)^^")/'-^ } • 

For the subsequent algorithmic computations the minimum polynomial MipO(Q((^„) of Cn 
over (Q will be needed since it shows how to replace certain higher powers of Cn by sums of 
lower ones. As it turns out in Proposition |S1 MipO(q(C„) is simply the following nth cyclotomic 
polynomial. 

Definition 1. The nth cyclotomic polynomial is given by 

Fn := lliX - C) , 

C 

where C runs over all primitive nth roots of unity in C 

Lemma 2. For n € N, one has: 

(a) Fn is monic and deg{Fn) = (f>{n). 

(b) UFd = X^-l. 

d\n 

(c) Fn^nx]. 

Proof See Chapter VI.3 of Lang (1993). □ 

Remark 4. Lemma 13 shows that we can compute the nth cyclotomic polynomial recursively 
by use of the Euclidean algorithm in 1i[X]. 

Proposition 3 (GauiJ). The minimum polynomial MipO(Q(C„) o/Cn over Q is the nth cyclo- 
tomic polynomial Fn. 
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Proof. By Definition Q Cn is a root of Now, note that MipO(Q(C„) is, by definition, the 
(uniquely determined) monic polynomial in of minimal degree having Cn as a root. Of 

course, it is a standard fact that deg(MipO(Q(^„)) = [K„ : (Q]; see Proposition 1.4 in Chapter 
V.l of Lang (1993). By Proposition^ one has [K„ : Q] = 4>{n), hence the result follows from 
Lemma 13 □ 

The final result of this preliminary section will provide a uniform finite upper bound on the 
number of 0„-equivalence classes in arbitrary grids for given X-ray directions. 

Proposition 4. If G is a torsion-free Abelian group of rank r, and H is a subgroup which is 
also of rank r, then the subgroup index [G : H] is finite and equals the absolute value of the 
determinant of the transition matrix A from any 'E-basis of G to any 'E-basis of H. 

Proof. See Chapter 2, Lemma 6.1.1 of Borevich & Shafarevich (1966). □ 

3. Model Sets 

Now we will first give a brief introduction to model sets and then we define the class of 
cyclotomic model sets that will be the underlying ground structure for the present paper. 

3.1. General Setting. By definition, model sets arise from so-called cut and project schemes. 
These are commutative diagrams of the following form; compare with Moody (2000) and see 
Baake et al. (2002) for a gentle introduction with many illustrations. 



(2) 



j^k ^ — U'' X H — > H 



U U lattice U dense 

_ 1-1 _ _ 
7r[L] < > L > 7Tint[L] 



Here, H is some locally compact Abelian group, vr and vTint are the canonical projections, and 
L is a lattice in R,*' x H, i.e., L is a discrete subgroup of R'^ x H such that the quotient group 

(R'' X H)/L 

is compact. Further, vTint [L] is a dense subset of H and the restriction of vr to L is assumed to be 
injective. Writing L := tt[L], one can define a map .*: L — > H hy x i — > iTi^t{Tr\2^ (x)) . Then, 
one has [L]* = 7ri„^[L]. If the map .* is injective, we denote the inverse of its co-restriction 
L — . [L]* by [L]* — > L. _ 

In the following we use the notation A°, A, dA for the standard topological operators 
interior, closure, and boundary of a set A in a locally compact Abelian group. 

Definition 2. (a) Given the cut and project scheme a subset W C H is called a 
window if / W° CWC W° and W° is compact, 
(b) Given any window W C H, and any t G R'^, we obtain a model set 

A{t,W) :=t + A{W) 

relative to the cut and project scheme by setting 

A{W) :={x e L\x* eW}. 
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Further, M.^ (resp., H) is called the physical (resp., internal) space and W is also 
referred to as the window of A{t,W). The map .*: L — > H, as defined above, is the 
so-called star map. 

For details about model sets and general background material see Moody (2000) and Baake 
& Moody (2000); see Baake et al. (2002) for detailed graphical illustrations of the projection 
method. 

Remark 5. The translation vector t in Definition |21 stresses an intrinsic character of model 
sets. While the structure model specifies the cut and project scheme k, H and L, and also 
the window W , a natural choice of the origin is usually not possible. 

Remark 6. Without loss of generality, we may assume that the stabilizer H\y of the window 
i.e., 

Hw ■■= {h(^ H\h + W = W}, 
is the trivial subgroup of H , i.e., Hw = {0}. Observe that the latter is always the case if H 
is some Euclidean space, i.e., if one has H = n'^ for some suitable d E N. Note further that 
the star map is a homomorphism of Abelian groups. 

The following remark collects some properties of model sets; for details see Moody (2000). 

Remark 7. In the following, for x G R'^ and r > 0, we denote by Br{x) the open ball of 
radius r about x. The model set A := A{t, W) C R'^ is a Delone set, meaning that A is both 
uniformly discrete (i.e., there is a radius r > such that every ball of the form Br{x), where 
X S R*^, contains at most one point of A) and relatively dense (i.e., there is a radius R > 
such that every ball of the form Br{x), where x G R'^, contains at least one point of A). Also, 
A has finite local complexity, i.e., A — A \s discrete and closed. (Note that A has finite local 
complexity iff for every r > there are, up to translation, only finitely many point sets (called 
patches of diameter r) of the form A n Br{x), where x £ R*^.) In fact, A is even a Meyer set 
(i.e., in addition. A — A is uniformly discrete). 

Further, A is aperiodic, i.e., has no translational symmetries iff the star map is injective. In 
fact, the kernel of the star map is the group of translational symmetries of A. 

If A is regular, i.e., the boundary dW of the window W has (Haar) measure in H, then A 
is pure point diffractive {cf. Schlottmann, 2000). If A is generic, i.e., [L]*n dW = 0, then A is 
repetitive. This means that, given any patch of radius r, there is a radius R such that any ball 
Br{x) in R'^ contains at least one translate of this patch; see Schlottmann (2000). If yl is both 
generic and regular, the frequency of repetition of finite patches is well defined, i.e., for every 
finite patch, the number of occurrences of translates of this patch per unit volume in the ball 
i?r(0) of radius r about the origin approaches a positive limit as r — > oo; cf. Schlottmann 
(1998). 

For the discrete tomography of aperiodic model sets, one additional difficulty, in comparison 
to the crystallographic case, stems from the fact that it is not sufficient to consider one pattern 
and its translates to define the setting. In particular, to define the analogue of a specific crystal, 
one has to add all infinite patterns that emerge as limits of sequences of translates defined in 
the local topology (LT). Here, two patterns are e-close if, after a translation by a distance of 
at most e, they agree on a ball of radius 1/e around the origin. If the starting pattern P is 
crystallographic, no new patterns are added; but if P is a generic aperiodic model set, one 
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ends up with uncountably many different patterns, even up to translations! Nevertheless, all 
of them are locally indistinguishable (LI). This means that eyery finite patch in A also appears 
in any of the other elements of the Ll-class and vice versa; see Baake (2002) for details. 

Remark 8. The entire Ll-class of a regular, generic model set yl(Ty) can be shown to consist 
of all sets t + A{t + W), with t G and r such that [L]* D d{T + W) = (i.e., r is in a 
generic position), and all patterns obtained as limits of sequences t + A{Tn + W), with all in 
a generic position; see Baake (2002). Each such limit is then a subset of some t + A{T + W), as 
r might not be in a generic position. In view of this complication, we must make sure that we 
reconstruct finite subsets of generic model sets. This will be reflected in Definitions El and d 
of Sectional 

3.2. Cyclotomic Model Sets. In the present paper we will study the discrete tomography 
of a special class of planar model sets, the cyclotomic model sets, which can be described in 
algebraic terms and have an Euclidean internal space. In the following let n G N \ {1, 2}. 

Before we formally introduce the cut and project scheme from which the cyclotomic model 
sets arise, let us consider some main ingredients. 

The elements of the Galois group G{¥in/Q.) (see Proposition come in pairs of complex 
conjugate automorphisms. Let the set {cti, • • • , o"(/,(ra)/2} arise from G(K„/(Q) by choosing 
exactly one automorphism from each such pair. Here, we always choose cxi as the identity 
rather than the complex conjugation. Every such choice induces a map 

: 0„^(r2)4^ 

through 

(Actually, .~ and the following map .* are defined on K„, but it is their restriction to On that 
is relevant here.) 

With the understanding that for cj){n) = 2 (i.e., n G {3,4,6}), the singleton 

(r2)^-i = (r2)° 

is the trivial (locally compact) Abelian group {0} each such choice induces a map 

.-o„^(ii^)^-S 

defined by .* := 0, if n G {3, 4, 6}, and 

z ! — > (cr2{z), . . .,a^{zj^ 

otherwise. Then, [On]~ is a Minkowski representation of the maximal order 0„ of K^, 
see Chapter 2, Section 3 of Borevich & Shafarevich (1966) and Theorem 2.6 of Washing- 
ton (1997). It follows that [Onf is a (full) lattice in IR^ x (]R2)9i{n)/2-i_ jj^^.^^ gj^^^g ^j^g ^^^^^ 
R2 X (i^2)0{n)/2-i -g Euchdean, this means that there are 4>{n) R-linearly independent vectors 
in IR? X (]R,2-)(/>(n)/2-i i;jg^Yj]^g j^j^g property that [On]~ is the Z-span of these vectors; compare 
Chapter 2, Sections 3 and 4 of Borevich & Shafarevich (1966). In fact, the set 

{r,(C„)~ ...,(Cn*(")-^r} 
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has this property; cf. Proposition [21 and Remark |31 Further, the image [On]* is dense in 
. This follows for instance from the existence of a Pisot number of (full) degree 
4>{n)/2 in see Chapter 2, Section 3 of Borevich & Shafarevich (1966) and Pleasants (2000). 
Multiplication by such a Pisot number in the physical space then translates via the map .* into 
a contraction in all directions of the internal space, as defined by the Q-span of the projected 
basis vectors of the lattice. 

Now, the cyclotomic model sets arise from cut and project schemes of the following form, 
where we follow Moody (2000), modified in the spirit of the algebraic setting of Pleasants 
(2000). 



r2 ^ ^2^^^2^^-l (r2)^-1 



(3) 



U 

1-1 

On ^ 



U lattice 
[On]- 



U dense 
[OnY 



As described above, one has 

[OnT = \ (z,{a2{z),. ■ ■ ,Cr4>{n) (z))) 

^ 2 



zeO, 



Recall that for n 7^ 3,4, 6 also the first inclusion On C involves a dense set. Now here 
is the definition of the class of cyclotomic model sets; for more details and related general 
algebraic settings, see Pleasants (2000). 

Definition 3. Given any window W C (IR,^)''^^")/^"^, and any i G R^, we obtain a planar 
model set 

An{t,W) ■.= t + An{W) 

relative to the above cut and project scheme ^ (i.e., relative to any choice of the set 
{cTj \ i £ {2, . . . , (p{n)/2}} as described above) by setting 

An{W) := {zGOn\z* £W}. 

We set 

ten^wc (R2)^-i is 1 
a window j 
Then, the class CAi of cyclotomic model sets is defined as 

CM:= [j M{On). 

neN\{l,2} 

Remark 9. The set A := An(t,W) C is aperiodic iff n ^ {3,4,6}, i.e., the translates of 
the square (resp., triangular) lattice are the only cyclotomic model sets having translational 
symmetries; compare Remark If, for a given n, A is both generic and regular, and, if the 
window W has m-fold cyclic symmetry with m a divisor of lcm(n, 2) and all in a suitable 
representation of the cyclic group Cm of order m, then A has m-fold cyclic symmetry in the 
sense of symmetries of Ll-classes. This means that A and the structure obtained by applying 



M{0„ 



Anit,W) 
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Figure 2. A central patch of the eightfold symmetric Ammann-Beenker tiling 
with vertex set Aab (left) and the /-image of vIab inside the octagonal window 
in the internal space (right), with relative scale as described in the text. 



an appropriate 'symmetry' are locally indistinguishable (LI); see Baake (2002) for details on 
the symmetry concept. 

3.2.1. Some Examples. We give five examples of cyclotomic model sets. The first two are 
periodic of the form yl„(0, W) G M{On) with n G {3, 4} (and hence W = {0}), while the last 
three are aperiodic cyclotomic model sets of the form An{0, W) G M.{On), with n G {5, 8, 12} 
(whence having an internal space of dimension 2). 

(a) The planar, generic, regular and periodic cyclotomic model set with 4-fold cyclic sym- 
metry associated with the well-known square tiling is the square lattice, which can be 
described in algebraic terms as tIsq := ^4(0, W^) = = O^; see Figure [21 

(b) The planar, generic, regular and periodic cyclotomic model set with 6-fold cyclic sym- 
metry associated with the well-known triangle tiling is the triangle lattice, which can 
be described in algebraic terms as vItri '■= ^3(0, W^) = O3; see FigureEl 

(c) The planar, generic and regular model set with 8-fold cyclic symmetry associated with 
the Ammann-Beenker tiling (Baake & Joseph, 1990; Ammann et ai, 1992; Gahler, 
1993) can be described in algebraic terms as 

ylAB := {zG08\z*gW}, 

where the star map .* is the Galois automorphism in ^(Kg/Q), defined by Cs ' — > Cfi 
and the window W is the regular octagon centred at the origin and of unit edge length, 
with orientation as in Figure |21 This construction also gives a tiling with squares and 
rhombi, both having edge length 1; see Figure |2I 
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Figure 3. A central patch of the tenfold symmetric Tubingen triangle tiling. 

If t G \{0} is chosen such that t + W is again in a generic position (this is true for 
almost alH G R^), the replacement of W hy t + W again leads to an Ammann-Beenker 
tiling. Moreover, the two tilings are locally indistinguishable (compare Remark |H1). 

(d) The planar and regular model set with 10-fold cyclic symmetry associated with the 
Tubingen triangle tiling (Baake et at, 1990a, b) can be described in algebraic terms 
as 

Aij^rj.^ := {z e 05\z* et + W} , 

where the star map .* is the Galois automorphism in ^(Ks/Q), defined by Cs ' — > Cf- 
Moreover, the window W is the regular decagon centred at the origin, with vertices in 
the directions that arise from the 10th roots of unity by a rotation through vr/lO, and of 
edge length t/\/t + 2, where r is the golden ratio, i.e., r = (\/5 + l)/2. Furthermore, 
t is an element of R^. Note that yll^rprp is not generic, while generic examples are 
obtained for almost all i G R^. Generic yl^prprp always give a triangle tiling with long 
(short) edges of lengths 1 (l/r, respectively); see Figure 01 Different generic choices 
of t result in locally indistinguishable (LI) Tiibingen triangle tilings (compare again 
Remark IHI). 

(e) The planar and regular model set with 12-fold cyclic symmetry associated with the 
shield tiling (Gahler, 1993) can be described in algebraic terms as 

4 := {z £Oi2\z* £t + W} , 

where the star map .* is the Galois automorphism in G(Ki2/Q), defined by C12 ' — > C121 
and the window W is the regular dodecagon centred at the origin, with vertices in the 
directions that arise from the 12th roots of unity by a rotation through vr/12, and of 
edge length 1. Again, t is an element of R^. Note that ylg is not generic, while yl| is 
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Figure 4. A central patch of the twelvefold symmetric shield tiling. 

generic for almost all t G R^. The shortest distance between points in a generic yl| is 
(\/3 — l)/-v/2. Joining such points by edges results in a shield tiling, i.e., a tiling with 
triangles, squares and so-called shields, all having edge length (\/3— 1)/\/2; see Figure^ 
for a generic example. Different generic choices of t result in locally indistinguishable 
shield tilings (compare again Remark [HI). 

4. Discrete Tomography of Planar Model Sets: Problems and Main Results 

4.1. Consistency, Reconstruction and Uniqueness. It is clear that each subset of the 
lattice is determined uniquely by one X-ray in an irrational direction. Therefore, the 
nontrivial classical problems of discrete tomography involve lattice directions, i.e. directions 
spanned by two lattice points. One now needs the correct analogue of lattice directions in the 
framework of cyclotomic model sets. 

Definition 4. Let n G N \ {1, 2}. 

(a) The elements of On \ {0} are called On- directions. 

(b) For an C^-direction o, we denote by Co the set of lines t + IRo with i G IR, while 
C'^" C Co is the set of module lines in direction o, i.e., the set of lines t + IRo in IR^ 
with t G On- 

Since every O^-direction is parallel to a non-zero element of the difference set 
An{t,W) — An{t,W) C On (Huck, 2006), the notion of On-directions is indeed the natural 
extension for cyclotomic model sets. 

Definition 5. Let n G N \ {1,2} and let F C IR^ be a finite set which lives on On, i.e., 
F <Z t + On, where t G R^. Furthermore, let o be an 0„-direction. Then, the {discrete 
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parallel) X-ray of F in direction o is the function 

XoF : ^ No := N U {0} , 

defined by 

XoF{£) :=card(Fnf). 

Remark 10. Obviously, XqF has finite support supp(XoF) (the set of hnes in direction o 
that pass through at least one point of F) and, moreover, 

XoF{l) = card(F) . 

^esupp(XoF) 

In view of the complications with limits indicated at the end of Section 13.11 we will make 
sure that we reconstruct finite subsets of generic model sets, i.e., subsets whose .*-image lies 
in the interior of the window. This restriction to the generic case is the proper analogue of 
the restriction to perfect lattices and their translates in the classical case. 

Definition 6. Let n G N\{1, 2}, let W C (]R,2)<^(n)/2-i ^ window (c/. Definition EJ, and let 
a star map .* be given, i.e., a map .*: On — > (]R^)'^(")/^~^, given by z i — > 0, if n € {3,4,6}, 
and given by z i — > {a2{z), . . . ,cr^(^n)/2{z)) otherwise (as described in Definition inj . Then, the 
elements of the subset 

{Anit,T + W°)\t Gie,T G (R^)^'^} 

of Ai{On), which are defined by use of the above star map are called W/XlCCO 

Remark 11. Let n G N\{1, 2}. Note that, if C (]R2)0(n)/2-i jg ^ window, then its interior 
W° is also a window. Note further that for n = 4 (resp., n £ {3, 6}) the set of W^XfCOn) *"Sets 
simply consists of all translates of the square lattice O4 (resp., triangular lattice O3). 

Definition 7 (Consistency, Reconstruction, and Uniqueness Problem). Let the data be given 
as in Definition El Further, let oi, . . . , Om be m > 2 pairwise non-parallel 0„-directions. The 
corresponding consistency, reconstruction and uniqueness problems are defined as follows. 

Consistency. 

Given functions ■ — > No, i G {!,••• whose supports are finite 
and satisfy supp(poi) C decide whether there is a finite set F which is 

contained in a W^XlCOn) ^'^^ satisfies Xq^F = po^, i G {1, . . . , m}. 
Reconstruction. 

Given functions poi ■ Co^ — > No, i G {!,..., m}, whose supports are finite 
and satisfy supp(poi) C Z^^", decide whether there exists a finite set F in a 
^M{0 ) ^^^^ satisfies Xq^F = po^, i £ {I, . . . ,m}, and, if so, construct 
one such F. 

Uniqueness. 

Given a finite subset F of a 1^X1(0 ) decide whether there is a different 

finite set F' that is also a subset of a W^^^ ^ and satisfies Xq.F = Xo^F' , 
i G {!,... ,m}. 
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Note that the parameter n, the directions Oj, and the window W are assumed to be fixed, 
i.e., are not part of the input. 

For results on the computational complexity of these problems in the lattice case (and 
the Turing machine as the model of computation), see Gritzmann (1997) and Gardner et al. 
(1999). 

4.2. The Decomposition Problem. Now we introduce the problem of how to decompose a 
grid (c/. Definition IHI) into translates of maximal Ori-subsets. Note that the crystallographic 
cases, namely, the triangular lattice and the square lattice, are included. 

Definition 8. Let n G N \ {1,2} and let oi, . . . , Om be m > 2 pairwise non-parallel On- 
directions. Moreover, let ■ — > No, i £ {I, . . . ,m}, be functions whose supports are 
finite and satisfy 

supp(poJ C . 

Then, the associated grid Gj^^ |jg|i^ is defined by 

m I 

G{p„jig{i,...,™}} := n u 

«=1 yGsupp(po.) 

Definition 9. Let n G N \ {1, 2}. We define an equivalence relation ~n on by setting 

X ^nV x-y £On. 

If X, y S satisfy x ~„ y, we say that x and y are equivalent modulo On- 

Definition 10 (Decomposition Problem). Let n G N \ {1,2}, and let oi, . . . , Om be m > 2 
pairwise non-parallel 0„-directions. The corresponding decomposition problem is defined as 
follows. 

Decomposition. 

Given functions ■ Cq^ — > Nq, i G {l,...,m}, whose supports are finite 
and satisfy supp(j)oi) C C^", compute the equivalence classes modulo 0„ in 
the associated grid G|p^_|jg|i 

Of course, this problem can be reduced to a polynomial number of membership tests in On- 
It is, however, not directly clear how these tests can be performed and, actually, the algebraic 
properties of On will be utilized. Also, later a uniform bound for the number of classes will 
be given that is independent of the X-ray data. 

Remark 12. The phenomenon of multiple equivalence classes modulo On in the grid occurs 
already in the classical lattice situation; see Figure [HI on the left. There, no translate of the 
finite subset of the square lattice (marked by the connecting lines) is contained in any of the 
other equivalence classes. Also, note the fact that exactly one of the three equivalence classes 
has 14 elements (the equivalence class marked by light grey), whereas the remaining two only 
have 13 elements; it follows that this equivalence class (which generates the same grid as the 
marked finite subset of the square lattice) would be the unique solution of the corresponding 
reconstruction problem associated with its point set. Hence, the problem of decomposing the 
grid into its equivalence classes modulo On is the first problem to be solved when dealing with 
the consistency or the reconstruction problem, also in the classical planar setting. 
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Figure 5. Grids arising from two Cn-directions: On the left, the grid gen- 
erated by the X-rays of a finite subset of a translate of O4 = 1? in the two 
non-parallel 04-directions (1,1) and (1,-2). The three equivalence classes 
modulo O4 are marked by different greyscales. On the right, the grid generated 
by the X-rays of a finite subset of a translate of tIab in the two non-parallel 
Cs-directions 1 and C| = i. The two equivalence classes are also shown. 

4.3. The Separation Problem. When dealing with the consistency, reconstruction and 
uniqueness problems defined above, it is clear from the definition of W^^^ ^ ^-sets that, given 
n E N\ {1,2, 3, 4, 6}, a finite set F of points in (u^)<t>(n)/2-i and a window W C (Il2)'^(n)/2-i^ 
we have to be able to decide whether F is contained in a translate of W°. This leads us to 
the following geometric separation problem for sets [F]* C t + W. 

Definition 11. Let d G N, let P,W C R'^, and let t G K'^. We set 

SwAP) ■.= Pn{t + w) 

and, further, 

SepwiP) ■■= {SwAP)\i^'^''} ■ 

Definition 12 (Separation Problem). 
Separation. 

Given a finite set P C R^, and a set C M!^, determine Sepi4/(P). 

Remark 13. Note that Sep^y(P) contains all subsets of P that are 'separable' from their 
complement (in P) by a translate of W. Trivially, one has p(^t + WiSt£p — W. It follows 
that 

(4) SwAP)={P^P\i^P-W}. 

We will frequently make use of the above equivalence, because it allows us to switch between 
a separable set Sw,tiP) and the set of translation vectors that makes it separable; see Fig. El 
for an illustration. 
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Figure 6. On the left: If we translate W by t, then {pi,P2} is a subset of 
t + W, but {ps} is not. On the right: The 'world of translation vectors'. The 
point t is contained in pi — W and p2 — W, but not in ps — W. Again, we see 
that S\Y,tiP) = {pi,P2}- (In this example, the window is centrally symmetric 
with respect to the origin, i.e., W = —W.) 



4.4. Main Algorithmic Results. In the following we apply the real RAM-model of compu- 
tation, see e.g. Preparata & Shamos (1985). Here each of the standard elementary operations 
on reals counts only with unit cost. 

Our first result shows that the decomposition problem can be solved efficiently. 

Theorem 1. The decomposition problem can be solved in polynomial time in the real RAM 
model. More precisely, it is of complexity 0{s'^), where s is the maximum of the cardinalities 
of the supports of the given X-ray data functions. 

The next result deals with the separation problem. 

Theorem 2. Let the window W be given as an intersection of finitely many halfspaces, i.e., 
W = {x \ Ax < b} with A G IR,'^'^ and 6 S E,'. {The parameters d. A, and b are not part 
of the input). Then, for any finite set P C W^, the problem of computing Sep\Yo{P) can be 
solved in 0(card(P)'^^"'^) operations. 

As a consequence of Theorems ^ and |2l we see that the standard tomographic algorithms 
that have been developed for the lattice case can also be extended to the tomography of 
cyclotomic model sets. 

Theorem 3. Let W be given as in Theorem\B Then the problems Consistency, Recon- 
struction and Uniqueness as defined in Definition\^ can be solved with polynomially many 
operations and polynomially many calls to an oracle that solves the same problem on subsets 
of the plane of cardinality O(s^), where s is again the maximum of the cardinalities of the 
supports of the given X-ray data functions. 

As a simple corollary we finally note that the case of two directions can be solved in 
polynomial time even for cyclotomic model sets. 

Corollary 2. When restricted to two On-directions and polytopal windows the problems Con- 
sistency, Reconstruction and Uniqueness as defined in Definition\^ can be solved in 
polynomial time in the real RAM-model. 
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5. Analysis of the Problems, Proofs and More Results 

In the following we give a detailed analysis of the problems introduced in the previous 
section, prove the assertions stated there and obtain more results on the way. 

5.1. Tractability of the Decomposition Problem. We will now show that the number 
of equivalence classes of a grid is uniformly bounded by a number that depends on the given 
directions but is independent of the X-ray data. This result will then allow us to prove 
Theorem ^ 

Definition 13. Let n S N\{1,2} and let 01,02 be two non-parallel On-directions. We define 
the complete grid 0^01,02} 

Gioi,o2} - n ( u ^) ■ 

Proposition 5. Let ?i € N\ {1, 2} and let oi, 02 be two non-parallel On- directions. Then, the 
complete grid G^^-^ g^} satisfies On C Gj^^ ^ o.i^'d ^^{01,02} -^{01,02}; where one sets 

and a, 13,^,6 £ On are determined by oi = a + /?Cn and 02 = 7 + 5Cn- 

Remark 14. Note that the linear independence of {01,02} and {l,Cn} over IR, implies that 
a6 — / 0. Also, by definition Mj^^ q^} is an On-module of rank 2 with basis 
{01/ {a6 — P^),02/{a6 — P^)}, and Proposition [3 shows that 

01,02 01,02 

Note further that there are examples where the inclusion G^o-^ g^j C M^^^ g^y is not an equality. 
This is due to the fact that M^^^ g^y depends on the scaling of oi and 02, while G{oi,02} does 
not. On the other hand, let 7 E On and consider the two non-parallel Cm-directions 1 and 
7 + Cn- Then, by (jSJ and Lemma Q one has 

On C G{i^^+^^} C M{i_^+(.^} = On 

and hence G^i^^_^_^^y = M|i ,^_,_^^} = On- Further observe that, for n G N \ {1,2,3,4,6}, the 
complete grid Gjo^ is a dense subset of the plane, because already its subset On has this 
property; cf. Remark ^ 

Proof of Proposition^^ The first inclusion is obvious by definition. 

Next, we claim that On C M^^^ g^y. Let o G On- By Lemmata), there are unique if,ip£On 
with o = If + ipCn- By the linear independence of {01,02} over R, there are unique x,y £ M. 
with xoi + yo2 = o. Hence 

{xa + y7 - v?) + {x(3 + y5 - V')Cn = 

and, using the linear independence of {1, Cn} over R, we get that xa+yy — ip = xl3+y6—tp = 0. 
In matrix notation, this means that 

a 1 \ ( ^ \ ^ ( f 
(3 6 J[y J W 
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Cramer's rule now implies that 

x = [^p5- V'7)/(a<5 - G On/{a5 - 13-/) 

and 

y = {aip- Pip)/{a6 - p-f) G On/{a6 - ^7) . 

This proves our claim. 

Finally, consider g G G'joj q^}- By definition, there are elements o',o" G 0„ such that 
{g} = {o'+ Roi ) n (o" + IR,02 ) . Moreover , there are unique x,y with g = o' +xoi = o" + y 02 . 
Hence, xoi + {—y)o2 = o" — o' £ On and, by the same calculation as above, we get that 
x,y £ On/ictS — f3'j). Together with our first claim, this shows that g £ -^^{01,02}- '-' 

Lemma 3. M^^^ ^^y is a Z-module of rank (pin). 

Proof. This is an immediate consequence of the facts that Mj^^ q^} is an o^-module of rank 2 
and On is a Z-module of rank (p{n)/2; see R emark ITU and Remark 01 □ 

The following lemma shows that Mjo^ q^}; and thus 0^^1,02}^ decomposes into finitely many 
equivalence classes whose number depends only on {01,02}. Note that the symbol U is used 
to indicate disjoint unions. 

Lemma 4. The subgroup index [M{o^ q^} • ^n] is finite. Hence, there are c G N and 
ti,t2, ■ ■ ■ ,tc G Mjoj such that 

where, without restriction, ti = 0. It follows that every subset G of Mjo^ ,,2} satisfies the 
decomposition 

G = [j. (Gn{ti + On)) . 

Proof. By Lemma 01 M^g^ g^y is a Z-module of rank (p{n). Moreover, M^g^ ^^y is torsion-free 
because it is a subset of the field K„; see Remark El But On is a Z-module of rank (p{n) as 
well; see Remark |S1 Now, Proposition 01 yields the result. □ 

Remark 15. By Proposition 01 the subgroup index 

[Miou02} On] 

equals the absolute value of the determinant of the transition matrix A from any Z-basis of 
-^{01,02} to any Z-basis of On- It follows that, given the Z-coordinates of oi and 02 with 

respect to the Z-basis {1, Cn, Cn; • • • > Cn^"^ ^} of On [cf. Remark (31), one is able to compute 
[M|o^ 02} ■ On\. Note that, for any ^ £On, one has 

[Mii,y+c„}:On] = l; 

see Remark fTH 

Remark 16. Let n G N \ {1,2} and let oi, . . . , 0^,^ be m > 2 pairwise non-parallel 0„- 
directions. For any instance of the corresponding decomposition problem, the associated grid 
G{po^\iG{i,...,m}} satisfies 

card(G{p^ji6{i^...,^}}) < 
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where 

s := max({card(supp(poJ) | i G {1, • • • ,m}}) . 
Since G|p^ C G{oi,02}j Proposition shows that the last part of Lemma 0] apphes 

In the following we assume that the elements of the supports of the Po^, i G {1, • • • 
are given in the form o + Roj for suitable o G 0„. Moreover, we assume that all o's and all 
Oi are given through their Z-coordinates with respect to the Z-basis {1, Cm Cm • • • i Cn^"^ ^} of 
On (cf. Remark 131). 

We now prove Theorem ^ which we restate in a rephrased form. 

Theorem 4. The decomposition problem can be solved with O(s^) many real number opera- 
tions. 

Proof. The algorithm performs the following steps. 

Step 1: By the proof of Lemmata), the Euclidean algorithm in the inductive 

computability of the nth cyclotomic polynomial F„ = MipoQ(Cn) {cf. Remark and Propo- 
sition 13), the proof of Proposition and the Gaussian elimination algorithm, we are able to 
compute the Q-coordinates of the elements of the grid G{po. with respect to 

the Q-basis 

Sn^ Sni ■ ■ ■ 1 Sn J 

of K„ {cf. Proposition Pi efficiently. 

Step 2: Since {1, Cn, Cm • • • ' Cn^"^ ^} is simultaneously a Q-basis of K„ and a Z-basis of On 
{cf. Proposition H and RemarklHJ, one has for all gO)?!) • • • )'Z0(n)-i G Q the equivalence 

^g-j 90 + 91 Cn H h g</,{n)-lCn^"^~"^ G On 

90,9l; • • • )90(n)-l G Z . 

By Step 1, the elements of Gip^. |ig{i,...,m}} are given in the form 

90 + giCn + • • • + 9<^(n)-lC^^"^"^ , 

where 90i 9ii • • • 1 9(/.(n)-i G Q. Now, proceed as follows: choose an arbitrary element g 
of Gjp^ and compute the Q-coordinates of the differences g — h with respect to 

{l,Cn,C^,...,Cn^"^"'}, where /i G |je{i,...,m}} \ {a}- By the above criterion (EJ, a fixed 
h lies in the same equivalence class modulo On as g iff all coordinates oi g — h are elements 
of Z. Iterate this procedure by successively removing the computed equivalence classes and 
proceeding with the remaining subset of the grid and an arbitrary element therein. 

We already saw in Remark [THl that the last part of Lemma applies to Gjp^ 
This immediately implies that Step 2 of this algorithm computes the equivalence classes of 
the grid modulo On in at most 

c:= [M{„^,„,} : 0„] G N 

iterations. The inequality card(G|p^ m}}) ^ •s^ {cf- Remark now completes the 

proof. □ 
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Remark 17. The proof of Theorem 0] indicates that we actually do not need the full strength 
of the real RAM-model of computation. Rather, a Turing machine model that is augmented 
for algebraic computations suffices, see e.g. Buchberger et al. (1982). Then, of course, the 
binary size of the input matters. 

5.2. Tractability of the Separation Problem. The problem Separation in its general 
form is interesting on its own and we show now how to deal with it for windows W that are 
open polyhedra, i.e., 

W = {x\ Ax<b} with A £ B.^""^ and 6 G , 

where d > 2 is a fixed constant. The ideas presented here can be generalized to semialgebraic 
sets, but we prefer to keep the exposition more elementary. Also, polytopal windows with 
A^-fold cyclic symmetry, where N is the function from (P), are most relevant for model sets. 
(Note that the windows underlying the examples in Section 13.2.11 are polytopes.) 

We will begin with some standard facts about hyperplane arrangements as they are needed 
to deal with Separation. See Edelsbrunner et al. (1986) for more information on hyperplane 
arrangements, and Agarwal & Sharir (2000) and Halperin (2004) for surveys that cover also 
more general classes of arrangements. 

Definition 14. For i G {1, . . . ,/}, let ai e \ {0}, /3j G R, and consider the sets Hi = 
{x I afx = (3i}. Then Hi is called hyperplane and 7i = {Hi, . . . , Hi} is a hyperplane ar- 
rangement in IR'^. The sign vector SV{x) of some point x G M!^ is defined component-wise 
via 

— 1 if afx < Pi 
SVi{x) := < if afx = /3i , l<i<l ■ 
^+1 ifafx>l3i 

For s G {±1,0}' with 

Cs := {x I SV{x) = s} ^ , 
Cs is called a (proper) cell of the arrangement Ti.. 

Remark 18. The cells of an arrangement are relatively open sets of various dimensions. In 
particular, a cell Cs with sign vector s is full- dimensional if and only if s G {±1}'- Of course, 
is the disjoint union of all the cells of a hyperplane arrangement. Figure gives some 
illustration. 

In view of their general relevance, hyperplane arrangements are well studied and also algo- 
rithmically well understood. 

Proposition 6. Let Ti = {Hi, . . . ,Hi} he a hyperplane arrangement in IR*^. There exists an 
algorithm that computes a set of points meeting each cell ofTi in 0(1'^) operations in the real 
RAM model. 

Proof. Cf. Theorem 3.3 of Edelsbrunner et al. (1986). See also Chapter 7 of Edelsbrunner 
(1987). □ 

The proof of Theorem |21 will now be based on the following observation that ties the sepa- 
ration problem to certain hyperplane arrangements. 
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Figure 7. A small example for a hyper plane- arrangement in the plane. The 
hyperplanes are given by Hi = {x \ ajx = Pi}, where oi = (—1, 1), (3i = —1, 
02 = (-1,4), P2 = 3, as = (0,1), (3s = h = (-2,-3), = -4. On 
the right, the cells are drawn schematically. The arrangement consists of six 
points, 16 one-dimensional cells (thick lines) and 11 full-dimensional cells (grey 
areas). Some sign vectors are given. Note that not all vectors in {±1,0}^ occur 
as sign vectors of cells; e.g., (0,0,0,0) and (-1,-1-1,-1,-1) are not realized. 

Lemma 5. Let P = {pi, . . . ,pq} be a finite set of points in M!^, let W = {x \ Ax < b} with 
A £ B}^'^, b e R', and let aj denote the ith row of A, 1 < i < I. For l<i<l and l<j<q, 
set 

H'f^ := {x I aJx = {Ap, - 6),} 

Further, set 

n{W, P) := {i^p^ I 1 < i < / , l<j<q}. 
Then, one has the following: 

(a) The set pj — W is an intersection of open half spaces defined by the h'^^ , • • • , H^^'' , more 
precisely 

Pj - W = {x \A^x > Apj - b} . 

(b) For each cell Cg of the hyperplane arrangement 7i{W,P) with sign vector s = {sij)ij 
the following implication is true: 

t,t'eCs =^ SwAP) = Sw,t'{P) ■ 
{Of course, the reverse implication is not true in general; see Figure\^) 

Proof. Part (a) follows from a simple computation: 

p — W = {p — x \ Ax < b} = {x \ A(p — x) < b} 

= {x \ Ax > Ap-b} . 

For (b), recall from (0]) that 

SwAP) = {pj\i<j<q , t£pj-w} 

for any t £ R^. Using (a) we see that t £ pj - W iS SVij{t) = ■■■ = SVij{t) = +1. Now, if 
t,t' £ Cg, we have SV(t) = SV{t') = s, concluding the proof. □ 

Here is a restatement of Theorem |2I 
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Figure 8. This shows in essence the same situation as in Fig. El on the right 
we added the arrangement Ti.{W,P). Here, Sw,t{P) = Sw,t'iP), but t and t' 
are in different fulldimensional cells of the arrangement Ti.{W,P), see Lemma 
El Therefore the inverse direction of the implication in Lemma (b) is not 
true. 



Theorem 5. Let W = {x \ Ax < b} with A G b!^'^ and b e B! . Moreover, let P = 
{pi, . . . ,pq} be a finite set of points in M!^. Then, Sep^(P) can be computed with the aid of 
at most 0{{lq)'^^^) operations in the real RAM model. 

Proof. Our algorithm to determine Sep^r{P) performs the following steps. 

Step 1: Compute {Apj — b)i for 1 < i < I and 1 < j < q, to specify the hyperplane arrangement 

H{W,P) from LemmajSl 

Step 2: Compute a set T of points meeting every cell of TLiyV^P). 
Step 3: For each of the points t G T obtained in 2, compute Sw,t{P)- 
Step 4-' Output the collection of all the Sw,t{P)- 

The correctness of this procedure follows directly from Lemma 

Now we show the complexity assertion. Step 1 needs no more than 0{lq) operations. Step 2 
requires 0{{lqY) operations by Proposition (HI For Step 3, we decide \i t ^ pj — W for each j. 
To this end we test if t satisfies the inequalities a^t > {Apj — b)i, l<i<l,l<j<q. This is 
done with 0{lq) operations. In total we do not need more than 0{lq + {IqYlq) = 0{{lqY~^^) 
operations. □ 

Remark 19. As the proof of Theorem shows, if the number of hyperplanes defining the 
window W is regarded constant, then 

card(Sepv^^o(P)) = O (^card(P)'^) . 

Remark 20. Theorem can be generalized to semialgebraic sets W . The corresponding 
algorithm is then based on an analogue of Proposition El in the semialgebraic world; see Basu 
et al. (1996) and Theorem 2 of Basu et al. (1997). 

5.3. On the Tractability of Consistency, Reconstruction and Uniqueness. As a con- 
sequence of Theorems ^ and |2l we can now prove Theorem (HI In the following we only deal 
with Consistency in detail; the proofs for the other two problems are similar. As TheoremOl 
states, we want to reduce Consistency to a problem in the classical (anchored) case. 
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The number m G N \ {1} of X-rays and the different directions oi, . . . , are of course 
fixed as usual. 

AnchoredConsistency. 

Given s € N and Po^ '■ — > No, i G {1, • • • with finite supports whose 
cardinahtes are bounded by s, and a finite set 5 C with at most points. 
Decide whether there is a set F contained in 5 which satisfies Xq^F = poi, 
i G {1,... ,m}. 

Now we show that for polytopal windows the problem Consistency for cyclotomic model 
sets can be reduced to AnchoredConsistency. Let A be an algorithm for solving An- 
choredConsistency. (In the following A acts as a black box subroutine for the reduction.) 

Theorem 6. Let W he given as in Theorem [B Then Consistency can he solved with 
polynomially many operations and polynomially many calls to A. 

Proof. The algorithm performs the following steps. 

Step 1: Check first the necessary condition that the cardinalities 

iesupp(poj 

coincide for each i. If this is the case, proceed with Step 2. Otherwise the instance is incon- 
sistent. 

Step 2: Compute the elements of the equivalence classes Gi of the associated grid 
G{p,^^...^p,^} modulo On, say 

■ c 

G{p.,,...,p.„} = Ui=i^^ ^ 
in terms of their Q-coordinates with respect to the Q-basis 

{i,c„,e...,cf")-n 

of K„ (c/. Proposition Pi . By Theorem this can be done efficiently. 

Step 3: For all i S {1, . . . ,c}, compute the .*-image [Gj]* of Gi. Note that we consider the 
star map here as a map 

This can be done efficiently. Due to the definition of W^^q^-^ ^-sets, a solution F C Gi for our 
instance must satisfy the condition 

(7) 3rG(R2)^-i : [FYct + W°. 

Recall that for n S {3, 4, 6}, condition Q is always satisfied and one can proceed with Step 4. 
Otherwise, compute the set Sep^yo([Gi]*). By Theorem [21 this can be done efficiently. Note 
that, for every i £ {l,...,c}, a subset F C Gi that satisfies condition Q has the property 
that [F]* C P for a suitable P G Sepp^/o ([Gj]*). Finally, compute, for all i £ {1, . . . ,c} and for 
all P G Sep^o([Gj]*), the pre-images S := [P]~* of P under the star map. This can be done 
efficiently. Note that, with the above restriction n ^ {3,4,6}, the star map is injective. 

Step 4-' If n £ {3,4,6}, consider the equivalence classes S := Gi, i G {!,..., c}, having 
the property that card(Gi) > A^. Otherwise, consider, for all i E {!,..., c} and for all 
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Figure 9. A (card(F)xcard(F))-grid (left hand side). The subsets of the grid 
that conform to the X-ray data correspond to permutation matrices, hence 
there are (card(-F))! of them. Assume that the internal space is also two- 
dimensional and that the star map .* acts as a bijection on the grid points. 
Then the grid in the physical space is mapped to a grid in the internal space 
(right hand side). (As concrete example of this situation, we recall the examples 
in Section [!l2.1l E.g., if n = 8 and the star map is defined by ^ then 
the grid in the physical space can be chosen as a finite patch of the Z-span 
linz;({l,^8}) (a lattice in K'^). It is in 1-1 correspondence to its image under 
the star map. If the window is chosen 'slim' enough, covering a 'diagonal' 
of the image grid, then there is only one of the (card(-F))! solutions (in the 
internal space) that can be covered by a translate of the window. Thus, there 
is also only one single solution in the physical space. 

P G Sepiyo([Gj]*), the subsets S := [P]^* of Gi having the property that card([S']"*) > N. 
Then apply A on each such S. The instance is consistent iff A reports consistency for one of 
the sets S. □ 

Note that for m = 2 a polynomial-time algorithm A is available; see e.g. Slump & Gerbrands 
(1982). There it is shown how to set up a capacitated network that admits a certain flow iff 
the consistency question has an affirmative answer. Points in the grid correspond to arcs in 
this network. If we want to forbid certain positions, we only have to cancel the corresponding 
arcs. Hence we obtain Corollary El for Consistency. 

The proofs for RECONSTRUCTION and Uniqueness are analogous. 

Remark 21. Note that the seemingly more natural approach to find subsets F C Gi first 
that conform to the X-rays, and check then whether is satisfied may lead to an exponential 
running time. In fact, FigurelHlgives a simple example with a unique solution but exponentially 
many subsets of the grid conforming to the X-ray data. 



It is a pleasure to thank U. Grimm for fruitful discussions. This work was supported by the 
German Research Council (Deutsche Forschungsgemeinschaft), MB and CH within the CRC 
701, and PC and KL through GRK 447/2. 



Acknowledgements 



26 



BAAKE, GRITZMANN, HUCK, LANGFELD, AND LORD 



References 

Agarwal, P. K. &; Sharir, M. (2000): Arrangements and their applications. In: Sack, J.-R. &; Urrutia, J. (eds.). 
Handbook of computational geometry. Amsterdam: North-Holland, 49-119. 

Alpers, A. & Gritzmann, P. (2006): On stability, error correction, and noise compensation in discrete tomog- 
raphy. SIAM J. Discrete Math. 20, 227-239. 

Alpers, A., Gritzmann, P. Ik, Thorens, L. (2001): Stability and instability in discrete tomography. In: 
Bertrand, G., Imiya, A. &; Klette, R. (eds.) (2001). Digital and image geometry. Lect. Notes Comput. Sci. 
2243, Springer, Berlin, 175-186. 

Ammann, R., Griinbaum, B. &; Shephard, G. C. (1992): Aperiodic tiles. Discrete Comput. Geom. 8, 1-25. 
Baake, M. (2002): A guide to mathematical quasicrystals. In: Suck, J.-B., Schreiber, M. &; Haussler, P. (eds.). 
Quasicrystals. An Introduction to Structure, Physical Properties, and Applications. Springer, Berlin , 17-48. 
http : //www . arxiv . org/math-ph/9901014 

Baake, M. &; Grimm, U. (2004): Bravais colourings of planar modules with iV-fold symmetry. Z. Kristallogr. 
219, 72-80. http://www.arxiv.org/math.C0/0301021 

Baake, M., Grimm, U. & Moody, R. V. (2002): What is aperiodic order? Preprint. 
http : //www . arxiv . org/math . HO/0203252 

Baake, M. & Joseph, D. (1990): Ideal and defective vertex configurations in the planar octagonal quasilattice. 
Phys. Rev. B 42, 8091- 8102. 

Baake, M., Kramer, P., Schlottmann, M. & Zeidler, D. (1990a): The triangle pattern - a new quasiperiodic 
tiling with fivefold symmetry. Mod. Phys. Lett. B4, 249-258. 

Baake, M., Kramer, P., Schlottmann, M. &; Zeidler, D. (1990b): Planar patterns with fivefold symmetry as 
sections of periodic structures in 4-spacc. Int. J. Mod. Phys. B4, 2217-2268. 

Baake, M. & Moody, R. V. (eds.) (2000). Directions in Mathematical Quasicrystals. CRM Monograph Series, 
vol. 13, AMS, Providence, RI. 

Basu, S., Pollack, R. & Roy, M.-F. (1996): On the combinatorial and algebraic complexity of quantifier 
elimination.. J. ACM 43, 1002-1045. 

Basu, S., Pollack, R. &; Roy, M.-F (1997): On computing a set of points meeting every cell defined by a family 
of poh'nomials on a variety. J. Complexity 13, 28-37. 

Borcvich, Z. I. & Shafarevich, I. R. (1966). Number Theory. Academic Press, New York. 
Buchberger, B., Collins, G.E., Loos, R. (eds.) (1982). Computing Supplementum. Springer, Wien. 
Cowley, J. M. (1995). Diffraction Physics. North-Holland, Amsterdam. 

Edelsbrunner, H. (1987). Algorithms in combinatorial geometry.. EATCS Monographs on Theoretical Com- 
puter Science, Vol. 10. Springer, Berlin. 

Edelsbrunner, H., O'Rourke, J. & Seidel, R. (1986). SIAM J. Comput. 15, 341-363. 
Fewster, P. F. (2003). X-ray Scattering from Semiconductors. 2nd ed. Imperial College Press, London. 
Gardner, R.J. & Gritzmann, P. (1997): Discrete tomography: determination of finite sets by X-rays. Trans. 
AMS 349, 2271-2295. 

Gardner, R. J., Gritzmann, P. k, Prangenberg, D. (1999): On the computational complexity of reconstructing 
lattice sets from their X-rays. Discrete Math. 202, 45-71. 

Gahler, F. (1993): Matching rules for quasicrystals: the composition-decomposition method. J. Non-Cryst. 
Solids 153-154, 160-164. 

Gritzmann, P. (1997): On the reconstruction of finite lattice sets from their X-rays. In: Ahronovitz, E. &; 
Fiorio, C. (eds.). Lecture Notes in Computer Science. Springer, London, 19-32. 

Gritzmann, P., Prangenberg, D., de Vrics, S. & Wiegelmann, M. (1998): Success and failure of certain 
reconstruction and uniqueness algorithms in discrete tomography. Int. J. Imaging Syst. and Technol. 9, 101- 
109. 

Gritzmann, P, dc Vrics, S. & Wiegelmann, M. (2000): Approximating binary images from discrete X-rays. 
SIAM J. Optim. 11, 522-546. 

Guinier, A. (1994). X-ray Diffraction in Crystals, Imperfect Crystals, and Amorphous Bodies. Dover Publi- 
cations, New York. 



DISCRETE TOMOGRAPHY 



27 



Halperin, D. (2004): Arrangements. In: Goodman, J. &; O'Rourke, J. (eds.). Handbook of discrete and com- 
putational geometry. 2nd ed., Discrete Mathematics and its Applications. Boca Raton, FL: Chapman &; 
Hall/CRC, 529-562. 

Herman, G. T. &; Kuba, A. (eds.) (1999). Discrete Tomography: Foundations, Algorithms, and Applications. 
Birkhauser, Boston. 

Huck, C. (2006): Uniqueness in discrete tomography of planar model sets. Preprint. 

Kisielowski, C., Schwander, P., Baumann, F. H., Seibt, M., Kim, Y. & Ourmazd, A. (1995): An approach 
to quantitative high-resolution transmission electron microscopy of crystalline materials. Ultramicroscopy 58, 
131-155. 

Lang, S. (1993). Algebra. 3rd cd., Addison- Wesley, Reading, MA. 

Moody, R. V. (2000): Model sets: a survey. In: Axel, F., Denoyer, F. & Gazeau, J.-P. (eds.). 
Prom Quasicrystals to More Complex Systems. EDP Sciences, Les Ulis, and Springer, Berlin, 145-166. 
http : //www . arxiv . org/math . MG/0002020 

Pleasants, P. A. B. (2000): Designer quasicrystals: cut-and-project sets with pre-assigned properties. In: 
Baake, M. &; Moody, R. V. (eds.). Directions in Mathematical Quasicrystals. CRM Monograph Series 13, 
AMS, Providence, RI, 95-141. 

Preparata, F. P. & Shamos, M. I. (1985). Computational Geometry: An Introduciton. Springer, New York. 
Schlottmann, M. (1998): Cut-and-project sets in locally compact Abelian groups. In: Patera, J. (ed.). Qua- 
sicrystals and Discrete Geometry. Fields Institute Monographs, vol. 10, AMS, Providence, RI, 247-264. 
Schlottmann, M. (2000): Generalized model sets and dynamical systems. In: Baake, M. &; Moody, R. V. (eds.). 
Directions in Mathematical Quasicrystals. CRM Monograph Series 13, AMS, Providence, RI, 143-159. 
Schwander, P., Kisielowski, C, Seibt, M., Baumann, F. H., Kim, Y. & Ourmazd, A. (1993): Mapping projected 
potential, interfacial roughness, and composition in general crystalline solids by quantitative transmission 
electron microscopy. Phys. Rev. Lett. 71, 4150-4153. 

Slump, C. H. &; Gerbrands, J. J. (1982): A network flow approach to reconstruction of the left ventricle from 
two projections. Comput. Graphics Image Process. 18, 18-36. 

Steurcr, W. (2004): Twenty years of structure research on quasicrystals. Part I. Pentagonal, octagonal, 

decagonal and dodecagonal quasicrystals. Z. Kristallogr. 219, 391-446. 

Washington, L. C. (1997). Introduction to Cyclotomic Fields. 2nd ed.. Springer, New York. 



Fakultat fur Mathematik, Universitat Bielefeld, Postfach 10 01 31, 33501 Bielefeld, Germany 
E-mail address: {mbaake,huck}@math.mii-bielefeld.de 
URL: http : //www. math. uni-bielef eld. de/ba£ike/ 

Zentrum Mathematik, TU Muxchex, Boltzmaxxstr. 3, 85747 Garching bei Munchen, Germany 
E-mail address: {gritzman,langfeld,lord}@ma.tuin.de 
URL: http : //www-m9 . ma . tum . de/ dm/homepages/ gritzmann/ 



